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W IPONiORING  MILITARY  ACTIVITY 

Office  of  Naval  Research 
Physics  Branch 
Washington,  D.  C.  20360 


It  was  shown  by  Newell  in  1962  that  the  extreme  value  and  first  passage 
time  distributions  of  various  types  of  common  Markov  processes  asymptotically 
approach  those  for  independent  random  variables.  In  view  of  the  great  simplifica- 
tion this  occasions  in  the  calculation  of  a number  of  important  properties  of 
Markov  processes,  it  is  clearly  of  interest  to  determine  in  some  detail  the  condi- 
tions on  both  the  time  and  space  variables  under  which  this  equivalence  holds,  in 
this  paper  we  investigate  and  establish  these  conditions  for  markov  processes 
described  by  the  Fokker-Pl anck  equation  and  express  them  in  simple  analytic  forms 
which  are  directly  related  to  the  coefficients  of  the  Fokker-Planck  equation.  To 
demonstrate  the  usefulness  of  these  conditions,  we  apply  them  to  two  represer^dti ve 
examples  of  Fokker-Planck  equations,  the  Ornstein-Uhlenbeck  process  and  the  Montroll- 
Shuler  model  for  harmonic  oscillator  dissociation.  It  is  shown  very  clear  y r 
these  examples  that  the  extreme  value  and  first  passage  time  distributions,  and 
thus  the  mean  extreme  and  mean  first  passage  times,  of  these  Markov  processes 
approach  very  closely  those  for  independent  random  variables  at  finite  values  of  the 
time  and  space  variables.  ’ 
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ABSTRACT 


It  was  shown  by  Newell  in  1962  that  the  extreme  value  and  first 
passage  time  distributions  of  various  types  of  common  Markov  processes 
asymptotically  approach  those  for  independent  random  variables.  In  view 
of  the  great  simplification  this  occasions  in  the  calculation  cf  a number 
of  important  properties  of  Markov  processes,  \ t is  clearly  of  interest  to 
determine  in  some  detail  the  conditions  on  both  the  time  and  space  variables 
under  which  this  equivalence  holds.  In  this  paper  we  investigate  and 
establish  these  conditions  for  Markov  processes  described  by  the  Fokker- 
Planck  equation  and  express  them  in  simple  analytic  forms  which  are  directly 
related  to  the  coefficients  of  the  Fokker-Planck  equation.  To  demonstrate 
the  usefulness  of  these  conditions,  we  apply  them  to  two  representative 
examples  of  Fokker-Planck  equations,  the  Ornstein-Uhlenbeck  process  and 
tne  Montroll -Shuler  model  tor  narrnomc  oscillator  association  It  is 
shown  very  clearly  in  these  examples  that  the  extreme  value  and  first 
passage  time  distributions,  and  thus  the  mean  extreme  and  mean  first  passage 
times,  of  these  Markov  processes  approach  very  closely  those  for  independent 
random  variables  at  finite  values  of  the  time  and  space  variables. 


Key  Words:  Markov  Processes,  Independent  Random  Variables,  Fokker-Planck 

Equation,  Ornstein-Uhlenbeck  Equation,  Extreme  Value  Distributions, 
First  Passage  Time  Distributions,  Mean  Maxima,  Mean  First  Passage 
Times . 
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1 . INTRODUCTION 

Most  previous  work  on  the  theory  of  extremes  in  probability  theory 
deal;  with  independent  random  variables  in  discrete  time.^  Physical 
processes,  however,  frequently  involve  dependent  random  variables  in  con- 
tinuous time.  It  is  thus  of  interest  and  importance  to  extend  the  theoiy 
of  extremes  to  dependent  variables  in  continuous  time.  This  paper  is  con- 
cerned with  such  an  extension  to  a specific  class  of  Markov  processes. 

The  theory  of  extremes  is  the  study  of  the  distribution  of  the  ex- 
treme values  (maximum  or  minimum)  of  a random  variable  within  a given  time 
interval.  The  distribution  of  extreme  values  is  closely  related  to  the 
distribution  of  first  passage  times  oc  the  random  variable  to  a prescribed 
boundary. 

Some  properties  of  extreme  value  distributions  for  Markov  processes 
have  been  investigated  previously . In  Ref.  3,  Newell  discusses  the 
asymptotic  extreme  value  distribution  for  one-dimensional  processes,  addressing 
himself  particularly  to  the  question  "as  to  which  of  the  common  types  of 
Markov  ptocesses  give  rise  to  extreme  value  distributions  like  those  obtained 
for  independent  identically  distributed  random  variables." 

The  calculation  of  extreme  value  distributions  and  their  moments 
for  Markov  dependent  variables  is  invariably  a much  more  difficult  task 
than  the  calculation  of  these  properties  for  independent  random  variables. 

It  is  reasonable  to  expect,  however,  that  the  extreme  values  of  a Markov 
dependent  random  variable  within  a given  time  interval  themselves  form  a 
set  of  independent  random  variables  if  the  time  interval  is  sufficiently 
long.  This  assumption  has  been  adopted  by  Gumbel^  and  other  workers. 

Thus,  for  instance,  when  considering  the  problem  of  maxima  of  levels  in 
a river  over  a period  of  many  years,  Gumbel  makes  the  assumption  that  the 
set  of  observations  of  these  maxima  in  a given  interval  of  time  forms  a 
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set  of  independent  random  variables  if  this  time  interval  is  sufficiently 
long,  say  one  year.  This  is  certainly  a reasonable  assumption.  Even 
though  the  maximum  heights  of  the  river  are  strongly  correlated  variables 
if  observations  are  made  within  ten  minute  intervals,  one  would  certainly 
expect  that  this  correlation  has  completely  disappeared  as  the  observation 
interval  is  extended  to  one  year.  It  is  of  interest  to  know  how  rapidly 
these  correlations  decay  in  a Markov  process  so  that  one  can  calculate, 
rather  than  guess,  time  intervals  for  which  the  assumption  of  independent 
random  variable  behavior  of  the  extreme  values  is  valid. 

A stochastic  process  with  independent  random  variables,  can  be 
characterized  by  the  simple  functional  form  of  the  probability  F(^,t) 
that  the  process  stays  below  5 throughout  the  time  t.  This  form  is 
F(c,t)  = exp  [-xU)t],  where  the  function  x(0  is  related  to  the  statistical 
distribution  of  tho  rando"1  variahies  and  is  hence  dependent  on  the  details 
of  the  process.  From  this  distribution,  moment  properties  can  be  calculated. 
The  mean  first  passage  time  to  5,  for  instance,  is  given  by  T^(0  = [x(0]-1. 
The  mean  maximum  of  the  process  is  also  simply  related  to  F(c,t).  The  pro- 
bability F(s,t)  and  its  moments  for  a Markov  process  are  much  more  complicated 
and  in  general  difficult  to  calculate.  After  a sufficiently  long  time  t 
and  for  sufficiently  large  values  of  the  variable  the  Markovian  probability 
function  F(^,t)  is  well  approximated  by  exp[ - ( ^) t] , the  simple  form 
characteristic  of  an  independent  random  variable  process.  The  function 
XoU)  depends  on  the  details  of  the  Markov  process  but  is  much  simpler  to 
evaluate  than  the  exact  probability  F(s,t).  Approximate  moment  properties 
can  again  readily  be  obtained;  for  instance,  the  mean  first  passage  time 

to  s is  T M)  * [x  U)]"1  • 

/ 1 0 
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It  is  thus  of  importance  to  determine  conditions  on  both  the  time 


r 

and  space  variables  for  which  extreme  value  distributions  and  their 
moments  in  Markov  processes  approach  the  behavior  characteristic  of  pro- 
cesses of  independent  random  variables,  since  calculations  are  then 

greatly  simplified.  In  this  paper  we  investigate  and  establish  such 

I 

conditions  for  the  class  of  Markov  processes  which  obey  a Fokker-Planck 
equation.  We  find  such  conditions  and  express  them  in  forms  which  can 
be  directly  related  to  the  coefficients  in  the  Fokker-Planck  equation. 

These  conditions  are  then  tested  for  two  particular  processes,  namely 
the  Ornstein-Uhlenbeck  (O.U.)  process^  and  the  Montrol  1 -Shuler  model  ^ 
for  harmonic  oscillator  dissociation  in  the  high  temperature  limit. 

It  is  shown  in  these  examples  that  the  extreme  value  and  first  passage 
time  distributions  of  these  Markov  processes  approach  very  closely  those 
for  independent  random  variables  at  finite  values  of  the  time  and  space 
variables. 

In  Section  2 we  present  the  necessary  definitions  and  general 
expressions  for  first  passage  time  and  extremgm  value  distributions  and 
their  moments.  In  Section  3 results  for  independent  random  variables 
are  briefly  reviewed.  In  Section  4 we  consider  the  general  Fokker-Planck 
equation  and  establish  conditions  for  approach  to  independent  random 
variable  behavior.  These  conditions  are  tested  on  the  Ornstein-Uhlenbeck 
equation,  for  which  detailed  numerical  calculations  are  presented  in 
Section  5,  and  on  the  harmonic  oscillator  dissociation  equation  in  Section  6. 
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2.  DEFINITIONS 

In  this  section  we  define  the  functions  and  processes  to  be 
studied  in  the  body  of  this  paper. 

We  limit  oui  considerations  to  one-dimensional  Markov  processes 
with  random  variables  X ( t ) defined  for  continuous  time.  Let  f(x,  t xQ)dx 
be  the  probability  that  X(t)  lies  within  (x,  x+dx)  without  ever  haviny 
crossed  x = c in  time  (C , t),  given  X(0)  = xQ.  As  will  be  seen  below, 
all  quantities  relevant  to  the  theory  of  extremes  can  be  expressed  in 
terms  of  the  probability  density  f(x,t|x  }. 


Define  a new  random  variable 


Z(t)  s max  | X(i),  0 * x s t|  . (II. 1) 

Then  the  cumulative  distribution  function,  defined  by 

I 

FU.t|xQ)  S ft  dx  f(x,t|x0)  (II. 2a) 

has  the  probabilistic  meaning 

FU,t|x0)  * Prob  |z(t)  < £|X(0)  = xQJ  (II. 2b) 

where  x = r denotes  a reflecting  boundary;  r may  be  finite  or  negative 
infinite,  with  r s xq  < £.  The  function  F is  useful  when  treating  maxima 
problems.  Minima  problems  can  be  treated  similarly  but  then  £ < xQ  s:  r 
with  r finite  or  positive  infinite. 

An  important  random  variable  is  the  time  t when  the  process  crosses 
x * £ for  the  first  time,  i.e.,  the  first  passage  time 

T(0  Emin  }t|X(t)  * £}  . (II. 3) 

Since  F( t 1 xq)  is  the  probability  that  X(t)  never  crosses  x = £ during 
the  time  0 £ t s t,  it  can  be  related  to  the  first  passage  time.  From 
Eqs.  (1 1. 2b)  and  (1 1. 3)  we  find 

FU,tjx0)  = Prob  |T(0  > t|X(0)  = xQ[  . (II. 4) 

Two  more  useful  definitions  are 

Y(£,t|x0)dS  e [1^  F(£,t|x0)  ]d£  (II. 5) 

s Prob  |£  < Z(t)  < £ + d£ I X (0)  = x | 

♦(£,t|x0)dt  = - [|y  F(£,t|x0)  ] dt  (II. 6) 

= Prob  1 1 < T (£)  as  t + dt  I X (0)  = xj  . 


The  functions  f,  F,  f,  and  * are  Green's  functions  for  processes 
with  arbitrary  initial  distributions.  They  can  be  averaged  over  any 
initial  distribution  q(xQ) : 

FU.t)  e J*J  dxQ  n(xQ)  F(j.,t|xo)  (II. 7) 

* ProbjT(^)  > tj  = Probjz(t)  < 

vU. t)  * I(r  dxQ  n(xQ)  t(4.t|x0)  = |^FU,t)  (II. 8) 

$U,t)  ■ /£  dxQ  n(x0)  i(t.t|x0)  F(t.t)  • (II. 9) 

Using  these  distribution  functions  one  can  now  write  down  moments, 
i.e.,  various  mean  properties,  such  as  for  instance  the  mean  maximum  of 
the  variable  X ( t ) in  a given  time  interval  (0,t)  or  the  mean  first 
passage  time  to  x = We  define  the  nth  moments  of  Z(t)  and  T(0  by 

in^lx0)  df-  f^.t|x0)  = rn  + n/r  dC  4n’1[l-F(4,t|xo)](n  J0) 

Tnuixo)  e J0  dt  tP  ♦u»tix0) = nJ0  dt  tP  1 F(^»tix0)  • (ii-11) 

Moments  with  respect  to  an  initial  distribution  n(xQ)  will  be  denoted  by 
Zp(t)  and  Tn(0*  The  quantity  Z^t)  is  then  the  mean  maximum  of  the 
process  in  the  time  interval  (0, t ) ; the  variance  of  the  maximum  in  this 
time  interval  is  Z2(t)  - Z^(t).  The  mean  first  passage  time  to  x = C is 
T]  (0  and  its  variance  is  T2(0  - T*U).  The  corresponding  conditional 

quantities  are  similarly  defined. 

We  will  also  consider  processes  for  which  f^x.t  x ) dx  is  the 
probability  that  the  random  variable  X(t)  lies  within  (x,  x + dx)  without 
ever  having  crossed  x = t £(£  i 0)  in  time  (0,t),  given  X(0)  = xQ.  This 
case  is  of  interest,  for  instance,  in  the  study  of  largest  deviations  of 
the  random  variable  from  equilibrium.  All  definitions  given  so  *ar  can 


r 


be  suitably  modified  to  cover  this  situation.  We  define  the  random 
variable 

Y(t)  e max  j | X (T ) j ,0  s T a tj  (11.12) 

as  the  greatest  excursion  of  the  random  variable  X from  zero  in  the  time 
interval  (0,t).  The  appropriate  cumulative  distribution  function  for  this 
problem  is 

F(c,t|xo)  = J dx  f(x,t|xo)  . (11.13) 

The  probabilistic  definitions  in  Eqs.  (II. 2b)  and  (II. 4)  hold  here  with 
Y(t)  as  defined  in  Eq.  (11.12)  and  T(0  given  by 

T(0  e min  | t | X (t)  = t cj  . (11.14) 

Equations  (II. 5)  - (II. 9)  and  (II. 11)  remain  unchanged  except  for  the 
replace"’?''*,  r - - \ in  thp  lnwer  limits  of  integration.  In  place  of 
Eq.  (1 1. 10)  we  now  have 

Yn(t| x ) E f »U.t|x  ) « n f dc  C ^ 

Jo  0 t (11.15) 

with  F(^,t|xQ)  defined  in  (11.13),  the  quantities  Y^(t|xQ)  and  Y^t) 

are  now  mean  extrema;  T^(5|xq)  and  T^(s)  are  mean  first  passage  times  to 

x * t £.  It  should  be  noted  that  in  the  frequently  occuring  situations 

in  which  F(^,t)  is  an  even  function  of  C,  the  moments  Tn(0  and  Yn(t) 

here  are  identical  with  the  moments  Tn(0  and  Zn(t)  discussed  earlier  with 

a reflecting  boundary  at  r = 0. 

We  restrict  our  discussion  to  Markov  processes  in  continuous  time. 
We  believe  that  many  of  our  conclusions  are  valid  for  a large  class  of 
master  equations  of  the  form 

|^-f(x,tixo)  = Jdx  'K(x,x  ' ) f ( x * , t j xq  ) (11.16) 


where  K(x,x')  is  the  transition  rate  from  state  x’  to  state  x. 
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However,  we  will  concern  ourselves  only  with  a special  case  of  Eq.  (11.16), 
namely,  the  Fokker-Planck  equation 


a_ 

at 


f(*,t!V  = It  i"  "V*)  f(x.t|xQ)  + 


(11.17; 


with  the  Initial  condition  f(x,0jxQ)  = 6 (x  - xQ)  and  appropriate 
boundary  conditions  to  be  discussed  in  Section  4. 


3.  RESULTS  FOR  INDEPENDENT  RANDOM  VARIABLES 


Since  it  is  one  of  the  objectives  of  this  paper  to  compare  the 
results  for  extreme  value  distributions  of  Markov  processes  with  those  for 
independent  random  variables,  we  collect  here  the  relevant  results  for 
the  latter  case.^'  ^ The  statistics  of  extremes  for  independent  random 
variables  are  traditionally  dealt  with  in  terms  of  variables  defined  for 
dis-ret1  times.  We  will  thus  state  results  in  terms  of  processes  in 
discrete  time  and  then  take  limits  to  continuous  time  where  appropriate. 

Let  , X£,  ...  , Xn  be  n sequential  observations  of  a process. 
The  Xj 's  are  taken  to  be  identically  and  continuously  distributed 
independent  random  variables.  Let  Z(n)  he  a new  random  variable  defined  by 


Z(n)  e max  J Xm , 1 £ m £ n]  (III  .1 ) 

in  analogy  with  Eq.  (II. 1).  In  further  analogy  with  the  equations  in 
Section  2 we  can  define  probabilities  and  probability  densities  such  as, 
for  example 

F ( C » n)  * Prob  |z(n)  < $J  (III .2) 

(cf . Eq-  (II -7)  and 

n)  dt.  = [|t-  FU,  n>]  d? 

•=  Prob  Js  < Z(n)  < s + dr|  (III  .3) 
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[cf.  Eq.  (II. 8)].  The  moments  can  then  be  determined  as  in  Section  2. 

Here  the  probabilities  and  probability  densities  are  independent  of  the 
initial  condition  since  we  are  dealing  with  independent  random  variables. 

Let  P(x)  be  the  cumulative  distribution  function  for  each  x-  and 
p(x)  be  its  probability  density  function: 


P(x)  h Prob(X  < x) 


(III. 4) 


P(x)  e P(x) 


(III. 5) 


The  function  P(x)  [or  p(x)]  is  all  that  is  needed  to  determine  every 
property  of  the  process.  Thus,  for  instance,  the  probability  F(c,n)  that 
the  maximum  of  n observations  is  below  £ is  equal  to  the  probability 


that  each  of  the  , i s 1,  ...,  n,  is  less  than  5 so  that 


FU.n)  ■ teu)]" 


(III. 6) 


from  which  it  follows  that 


*U.r)  = n[P(0]P  pU) 


(III. 7) 


Me  ents  can  alc~  ie  given  immediately.  L:t  At  be  the  time  between  obser- 
vations and  let  t - iAt  be  the  time  of  the  nth  observation.  Time  is 
introduced  here  so  that  comparison  with  continuous  time  processes  can  be 
m?de  later;  we  could  of  course  continue  our  discussion  in  terms  of 
"• bservation"  or  "step  number"  and  never  introduce  time.  The  mean 
first  passage  time  to  x =£[cf.  Eq.  (II. 11)]  is  now  given  by  the  infinite 


series 


^(e)  = At | [l  - PU)]  +’ 2 PU)  [l  - PU)] 
+ 3 P2U)  [l  - PU)3  + ...| 


At 

T - p(0 


(in. 8) 


For  the  independent  process  we  define  an  upcrossing  rate  xU) 
so  that  xUUt  is  the  probability  that  the  process  crosses  the  level  x = s 
from  below  in  time  At.  Then 

x(f,)  At  = 1 - P(d  (III. 91 

and  the  probability  FU.t)  that  the  process  stays  below  5 throughout 
the  time  t = nAt  can,  instead  of  Eq.  (III. 6), be  written  as 

FU.t)  - [1  - xU)  at]"  . (111.10) 

In  the  limit  of  many  observations  such  that  n-*-®,  At-*-0,  nAt  s t, 

F(e,t)  — e'x(5)t  (HI- ID 

In  this  limit,  the  first  passage  time  moments  are  then  given  by  [cf . Eq. 

(11.11)] 

x U) 

and  the  mean  maxima  are  [cf.  Eq.  (II. 10)] 

Zn(t)  = rn  + n £ d?  cn‘1[l-e’X(0t]  . (III. 13) 

The  results  of  Eqs.  (III. 11)  and  (III. 12),  summarized  as 

F(ttt)_^e't/TlU)  (III. 14) 


are  characteristic  of  independent  processes.  The  form  of  *U)  or  T^U) 
of  course  depends  on  the  particular  process  one  is  considering,  i.e., 
on  PU). 

Equations  (III. 11),  (III. 12),  and  (III. 14)  also  hold  with  some 
changes  in  definitions  when  one  considers  first  arrival  at  x = -4.  In 
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I 


place  of  Eq.  (III. 4)  one  now  has 


P ( x ) ■ Prob 1 < x! 


(III. 15) 


|i  n\\  ' 

The  quantity  x(c)  At  is  the  probability  that  the  process  crosses  the 
level  x * £ from  below,  or  the  level  x * from  above,  in  time  At.  The 


moments  Yp(t)  now  are 


Y (t)  ■ n f dr 
•Jo 


(III. 16) 


Depending  on  the  form  of  the  cumulative  distribution  function 
of  Eq.  (III. 4),  and  hence  of  the  mean  first  passage  time  T^U).  one  obtains 
In  the  limit  as  t-*-“ (equivalent  to  n-*-“)  one  of  the  three  familiar 
stable  asymptotic  distributions  of  F(s,t)^*^.  The  most  familiar  of 

(13) 

the  three,  known  as  the  asymptotic  distribution  of  the  exponential  class, 
is 


with 


lim  F(f,  ,t)-e‘ 

t— ► °o 


Z *“n  - 6n> 


(III. 17) 


(III. 18) 


and  a ,6  defined  by  the  relations 
n n 


P(en)  = 1 • b an  “ " 


(III. 19) 


x = 6. 


This  asymptotic  distribution  holds  for  a wide  class  of  probability 
density  functions  p(x)  for  which,  as  x-*-« 


P"(x)  P*  (x)  P""(x) 


(III. 20) 


■11 


r 


Exponential  and  Gaussian  density  functions  are  members  of  this  class. 

The  quantity  en  is  the  expected  extreme  value  o X and  becomes  equal  to 
the  most  probable  extreme  value  when  n *> . The  quantity  an  is  a measure 
of  the  variation  of  3n  with  n.  The  mean  first  passage  time  to  x = 4 
corresponding  to  the  asymptotic  distribution  (III. 17)  is 


Vs) 

At 


1 -e 


-e 


-Z 


•1 


— eL  for  large  Z 
,(3,4) 


(III. 21) 


It  has  been  shownv'',',/  that  in  the  limit  c—  “,t— ® the  distribution 
F(s,t)  for  Markov  processes  tends  to  the  asymptotic  independent  random 
variable  form,  Eqs.  (III. 11)  and  (III. 12).  We  will  show  in  this  paper 
by  detailed  analysis  of  some  specific  Markov  processes  that  first  passage 
time  and  extreme  value  distributions  and  moments  are  given  to  a very  good 
approximation  by  the  indep  ndent  variable  formulae  also  for  finite  times 
ana  finite  boundaries.  The  orileria  for  this  ayi  ccnient  anu  the  ci'i'ui  3 
involved  will  be  discussed  in  detail  in  the  next  few  sections. 

4.  SOLUTIONS  OF  THE  FOKKER-PLANCK  EQUATION  WITH  BOUNDARIES 

In  this  section  we  formulate  the  solution  of  the  Fokker-Pl anck 

equation 

f (x ,t | xQ)  = 1^-  J-m^xJf  (x,t|x0)+£  ^ [m2(x)  f (x, t|  xq)]  J 

(IV. 1) 

in  the  presence  of  boundaries.  In  Eq.  (IV. 1)  the  coefficients  m^(x)  and 
m2(x)  are  given  by  the  conditional  expectation  values  of  the  change  aX 
of  the  random  variable  X in  time  At: 


E[aX | X ( t ) = x]  = m]  (x)  At  + o (a t ) 

E[ (AX )2 1 X( t ) = x]  = m2(x)  At  + o ( A t ) 
E[(AX)r|X(t)  = x]  = o(At),  r i 2 


(IV. 2) 
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The  probability  density  also  satisfies  the  backward  Kolmogorov  equation 


It  f<x't|*o> 

■ mi(x 0>  sr  f(x-tlxo)  * ?m2(xo)  r?f(x'tlxo)  • (lv-3) 

c xo 

There  exists  a large  literature  dealing  with  equations  of  the 
form  of  Eq.  (IV.  1).  Most  of  the  literature  can  roughly  be  divided  into 
three  (overlapping)  categories f ^ ^ One  category^ deals  with  the 
Ornstein-Uhlenbeck  (O.U.)  process  [m^(x)  = -px,  m^x)  = 2D,  p and  D 
constant]  as  a description  of  Brownian  motion,  where  the  velocity  distri- 
bution of  Brownian  particles  is  the  distribution  of  interest!  The  O.U. 

equation  is  solved,  but  not  for  extreme  value  distributions,  and  only 

(2-6  1 7-22 ) 

specialized  boundary  conditions  are  considered.  The  second  category'  " ’ ' 

deals  with  first  passage  time  or  extreme  value  distributions  of  diffusion 
processes.  In  this  group  are  uie  paper  which  are  Couccrncu  wiih  formal 
rather  than  detailed  solutions  of  Eq.  (IV. 1)  with  general  boundary  conditions 
and  with  asymptotic  properties  of  these  distributions.  Finally,  there 
are  many  papers  concerned  with  specific  applications  that  involve  various 
particular  forms  of  Eq.  (IV.l),  for  which  we  only  give  a few  typical  refer- 
ences. (9.23-25)  Qf  these  many  papers  only  the  work  on  the  Fokker- 
Planck  equation  described  in  refs.  3,  4 and  5,  deals  with  the  problem  of 
the  asymptotic  equivalence  of  the  extremum  properties  of  the  F-P  equation 
and  independent  random  variable  processes.  The  recent  numerical  work  of 
Keilson  and  Hiss^26^  on  passage  time  distributions  for  the  O.U.  process  is 
also  related  to  the  calculations  presented  here. 

k numL  -r  of  different  mathematical  techniques  have  been  used  to 
solve  Eq.  vlv.l)  and  to  find  properties  related  to  f(x,t|xQ).  The  two  most 

1 We  shew  in  Appendix  A that  the  O.U.  process  is  equivalent  to  a stationary 
auto-regressive  process  of  first  order(in  the  terminology  of  time  series 
analysis)  under  the  assumption  of  continuous  time  and  delta  correlated 
Gaussian  noise. 
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common  ones,  are  Laplace  transforms^18, 19,26^  and  eigenfunction  expansion.'3,20^ 
The  former  is  inconvenient  for  numerical  work  because  of  difficult  inverse 
transforms  that  must  be  performed,  and  becomes  tractable  only  for  certain 
specialized  boundary  conditions.  We  therefore  chose  the  method  of  eigen- 
function expansions. 

Equation  (IV. 1)  is  a parabolic  partial  differential  equation  whose 
"olution  can  be  expressed  as  the  eigenfunction  expansion 


V lMx)  ^Mxo^  -x  t 

f ( x , 1 1 x ) = 2-  p(x)  e k 

0 k=0  Nk 


(IV. 4) 


Here  the  weight  function  p(x)  is  the  equilibrium  solution  of  the  Fokker- 
Planck  equation  (IV. 1)  with  or  without  reflecting  boundaries,  but  without 


absorbing  boundaries, 


, . lx  2m,(x') 

p(x)  = — Mi  exp  \ dx'  i c i~r~  (IV. 5) 

with  h ( a ) determined  by  the  normalization  condition 

ao 

^ p (x)  dx  = 1 • ( IV .6) 

a 

The  constant  a is  left  undefined  at  this  point,  and  will  be  chosen  later 
to  conform  to  the  particular  boundary  conditions  to  be  considered.  The 
eigenfunctions  Uk  and  eigenvalues  *k  satisfy  the  differential  equation 


(IV. 6) 


i d\(x)  dUk(x) 

1 m2(x)  ■ — ^ +ml(x)  ~ + Xk  Uk(x)  = 0 


(IV. 7) 


In  addition  two  boundary  conditions  for  Uk(x)  must  be  specified.  The 
normalization  constants  Nk  also  depend  on  the  boundary  conditions. 

To  deduce  some  properties  of  the  eigenvalues  *k,  it  is  convenient  to 
rewrite  Eq.  (IV. 7)  in  the  self-adjoint  form 
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where 


dUk(x) 


+ Ak 


(IV. 8) 


e 2 m2(x)  p(x) 


(IV. 9) 


Since  this  equation  is  of  the  Sturm-L iouvi 1 le  type,  the  eigenvalues  *k 
are  real  and  non-degenerate.  For  m2(x)  > 0,  which  is  required  in  order 
that  the  F-P  equation  (IV. 1)  describe  a physical  system,  the  roots  are 


non-negative  and 


X < X,  < Xo  < • • 
0 I c 


(IV. 10) 


a)  Reflecting  boundary  at  x = r,  absorbing  boundary  at  x * £. 


We  impose  a reflecting  boundary  condition  at  x = r,  where  r may 
be  finite  or  negative  infinite.  For  definiteness  we  consider  here  maxima 
problems;  for  minima  problems  r is  finite  or  positive  intimte.  The 
boundary  condition  is 


L pn2(x)  ’) 

|ST  p-  f<x-tl>'o)J-ml(x)f(x-tlxo>Jx.r  ■ 0 


(IV. 11) 


for  all  t - G.  This  condition  insures  (provided  m^  and  m2  are  physically 
reasonable)  .'.at  the  process  is  conservative,  i.e.,  for  t ® 


$ Ck  c(x,t  x J = 1 


(IV. 12) 


The  second  boundary  condition  is  obtained  directly  from  the  definition 
of  f ( x , t | xQ ) dx  given  ir  Section  2 as  the  probability  that  the  random 
variable  X(t)  lies  within  (x,  x + dx)  without  ever  having  crossed  x = s 
in  time  (0,t),  given  X'O)  = xQ.  We  can  thus  impose  the  absorbing  boundary 


condition 


f ( x , 1 1 x ) I = 0 

'x=c 


(IV. 13) 


I 


r 


f 


It  should  be  stressed  here  that  the  condition  at  x . t does  not  necessarily 

It  sfmp1y  represer;$  the  #aiue 

{ °f  X “ “hiCh  “e  "watch”  f°r  the  first  passage,  or  crossing,  of  the 
random  variable  X(t).  The  initial  value  *p  of  X(t)  is  bounded  by  r * *Q  < 

-c.  tne:e  cornua  .■  .o-duions  we  take  a - r as  the  lower  limit 

for  the  no-Mli.ac,  r equal, (!».$,.  The  normalization  constants  N. 
in  Eq.  (IV/.71)  are  erven  by 


6k,mNk  = dx  °vx)  Uk(x)  Um(x) 

The  boundary  conditions  on  the  eigenfunctions  are 


dU.(x) 


e(0  U,  U)  = o 


(IV. 14) 


(IV. 15) 


(IV. 16) 


b ) Absorbing  boundaries  at  x - 

We  now  consider  the  case  of  symmetric  (about  ,.0)  absorbing  bound, 
anes  at  x = *4.  The  boundary  conditions  are. 


T(x,t|x 


'x  * 


(IV. 17) 


and  the  normal ization  constants  N,  are 


6k,m  Nk 


dx  p(x)  Uk(x)  U (x) 


(IV. 18) 


a • - - in  Eqs.  (IV. 5)  and  (IV. 6).  The  bounda 


ry  conditions  on  the 


e'icer  unctions  are 


o(*0  u,  (to  = 0 


(IV. 19) 


* 


c)  First  Passage  Times. 


To  obtain  first  passage  times  for  the  boundary  conditions 
considered  above  we  utilize  the  backward  Kolmogorov  equation  (IV. 3).  The 
first  passage  time  moments  Tn ( C I xq ) satisfy  the  differential  difference 
equations,  obtained  from  fcq.  (iv.3), 


V5|xo>  * 1 


<1  T„(C|«J 
dx 


MO  " 0 


IT  (C|x.) 


?"'2'V  tml<xo>  "d»  = -nTn-l(?ixo> 


(n  = 1,2,...)  (IV. 20) 


with  appropriate  boundary  conditions'27^  For  first 


passage  to  x = e with 


a reflecting  boundary  at  x = r,  the  boundary  conditions  on  T are 

n 


V^o> 


= 0 


x ■£ 
o 


(IV. 21) 


ajrTn(5lxo)|  =0 

o 1 x =r 

o 

For  first  passage  to  x = ±s,2 


(IV. 22) 


T"({|xo,L  5 0 


V*  < 


(IV. 23) 


Rewriting  Eq.  (IV. 20)  in  self-adjoint  form  and  integrating  twice  yields 

S z 

rn^|Xo^  = jix  ffzy  | n Tn-i ^ ly ) (iv.24) 


with  c = r for  boundaries  at  x = r , £(case  4a),and  c = 0 for  boundaries 


In  the  interest  of  simplicity,  when  dealing  with  two  absorbing  boundaries 
we  will  henceforth  assume  that  m^(x)  is  odd  and  ny/)  is  even. 
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at  x = -c  (case  4b). 

To  obtain  the  mean  first  passage  time  T^(0  for  the  initial  distri- 
bution rj(xQ)  = p(x  = xQ)  we  multiply  Eq.  (IV. 24),  with  n=  1,  by  p(xq)  and 
integrate  between  the  boundaries.  With  boundaries  at  r and  S, 


where 


(IV. 25) 


(IV. 26) 


With  boundaries  at  *£, 

T,(0  - ) <tz  (IV- 2?) 

Equation  (IV. 25)  reduces  to  tq.  (iv.^/j  when  r = u,  as  noted  earlier, 
d)  Asymptotic  Results. 

We  now  present  asymptotic  results  for  the  cumulative  distribution 
functions  FU,t|x  ) and  FU.t),  the  first  passage  time  moments  Tn ( C) , and 
the  extremum  moments  Zn(t)  and  Y (t)  for  the  Fokker-Planck  equation  dis- 
cussed above.  It  is  useful  to  define  new  functions  A^)  and  Bk(^)  by 

ao 

. » ) ■ Z AkldlUxJeV  (IV. 28) 

0 k=0  K K 0 


F(t.t)  = X C.  (c  )e‘  Ak.1 


(IV. 29) 


k=0 

from  which  it  follows  that  [see  Eqs.  (II. 2a),  (11.13)  and  (IV. 4)] 


Ak(c)  = n^hr  J> 


dx  o(x)  Uk(x), 


(IV. 30) 


with  b = r or  depending  upon  the  nature  of  the  boundaries. 
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(IV. 31) 


We  consider  in  this  paper  only  initial  probability  densities  n(x)  whi 
are  identical  with  the  equilibrium  weight  function  p(y)  defined  in  Eq. 

(IV. 5).  It  then  follows  that 

t 2 

BkU)  = AkU)  J dx  Vx^  = Vf^Ak^ 

i)  Cumulative  Distribution  Function^ 

For  large  r,  a U)— 0.  With  the  Pr°Per  choice  of  normalization, 
U (xo)-l  for  xQ«U  andN0(r)-l.  In  these  limits,  indicated  by  the 
symbol  ~,  with  one  reflecting  and  one  absorbing  boundary  it  follows  from 
Eqs . (IV. 30)  and  (IV. 31)  that 

oo 

A dx  p(x)  * 1-j  dx  p(x) 

0 L c 


and 


B0(0- 


I dx  p(x) 
LJr 


'2 

® 12 
c 

r 

1 - J dx  o(x) 

- 

C 

It  X - 

: one  finds 

(IV. 32) 


(IV. 33) 


C f 

A U)~  f dx  p(x)  = 1-  2 J dx  p(x) 

o f 


and 


Bq(0- 


ro 

oo 

f . 

1 dx  p(x) 

= 

1 - 2 J dx  p(x) 

J 

4 

(IV. 34) 


(IV. 35) 


It  now  follows  from  Eqs.  (IV.28)  and  (IV.29),  for  large  5 and  t,  independently  of 
boundary  conditions  and  initial  conditions,  that 

F((.tl«o)-F(l.O~«-*oU,t  (1V'36) 

• , (3-6) 

a result  which  has  been  derived  previously. 
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Since  8^(0  - 0 for  all  k and  (,  [cf.  Eq.  (IV. 31)]  while  Ak(OUk(x0) 
may  be  positive  or  negative,  one  can  be  more  specific  about  the  approach 
to  asymptotic  behavior  of  F(s,t)  than  one  can  for  F(4,t|xQ).  From  the 
initial  condition 

F U »t  = 0)  = X B,  U)  s f dx  p (x)  = e(0  (IV. 37) 

k=0  K Jb 

and  Eq.  (IV. 10)  and  Eq.  (IV. 29)  it  follows  that  for  all  (.  and  t 

B (c)e"Xo^^t  £ F (c , t>  sc(c)  e'Xo^^  . (IV. 38) 


Hence  Eqs.  (IV. 33),  (IV. 35)  and  (IV. 37)  imply  that,  for  all  times  t, 


FU.t)  « e~\>U)t 
for  £ sufficiently  large  so  that 

OD 

J p ( x ) dx  « 1 

it  should  be  noted  that  the  condition (IV .39a)  impl ies 


(IV. 39) 


(IV. 39a) 


B0U)  - eU)  - 1 


(IV. 39b) 

from  v>  * ch  (IV. 39)  follows  immediately. 

The  relations  (IV. 36)  and  (IV. 39)  are  identical  in  form  with 
F(c,t)  independent  random  variables  as  given  ,n  Eq.  (III. 11).  It  is 
in  .nis  sense  that  the  extremum  properties  of  Markov  processes  approach 
those  of  independent  processes,  even  for  finite  values  of  the  variables. 
For  the  Markov  process,  as  is  the  case  for  independent  random  variables, 
the  specific  form  of  ^ (£)  depends  on  the  details  of  the  process  under 
consideration. 
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i i ) First  Passage  Time  Moments 

F.om  Eqs . (Ii.11)  and  (IV. 29)  it  follows  that  the  first  passage 
time  moments  are  given  by 


k=0 


Bkt» 

*!<«> 


(IV. 40) 


Asymptotically,  when  Fqs.  (I''. 39)  and  (IV. 39a)  are  applicable,  this  reduces  to 


T (4)  — - — (IV. 41) 

" *"(4) 

which  is  identical  to  the  result  of  F'.  (II 1. 12)  for  independent  random 
variables.  In  particular,  Aq(0  is  again  the  reciprocal  of  the  mean  first 
passage  time  to  £ or  *£,depending  upon  the  boundary  conditions.  It  should 
be  noted  that  for  fixed  £,  Eq.  (IV. 41)  is  more  accurate  as  the  order  n 
of  the  moment  increases  since  the  rat  0 of  the  (k  + 1)^  to  the  k^  term 
in  the  sum  (IV. 40)  decreases  with  increasing  n.  inis  is  a consequence 
of  the  fact  that  the  accuracy  of  the  approximation  (IV. 39)  improves  for 
fixed  £ as  t increases.  This  is  an  important  observation  for  both  Markov 
and  independent  random  variables  since  the  mean  first  passage  time  by 
itself  is  a very  inaccurate  measure  of  the  time  dependence  of  any  process 
with  an  exponential  cumulative  distribution  of  the  form  (IV. 39). 

Indeed,  in  the  asymptotic  limits  where  (IV. 39)  and  hence  (IV. 41) 
are  valid,  the  dispersion  of  the  mean  first  passage  time  is  unity,  i.e. 


1 iin 
£ - * 


[( 


T2(£)  - T] ( c) 


f 


VO' 


(IV. 41a) 


The  utility  of  these  results  [Eq.  (IV. 40)  et.  seq.]  lies  in  the 
simplicity  with  which  one  can  calculate  Aq(0  = 1/^(0  from  Eq.  (IV. 25) 
or  Eq.  (IV. 27).  This  in  turn  permits  one  to  obtain  very  *-eadily  the 
asymptotic  distribution  function  (IV. 36)  and  (IV. 39). 


1 

I 
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i i i ) Extremum  Moments 

It  is  more  difficult  to  make  general  statements  about  the 
asymptotic  behavior  of  extremum  moments  than  of  first  passage  time  moments. 
The  latter  involve  integrations  over  time  of  a function  with  a simple  time 
dependence  which  is  known  for  all  t.  On  the  other  hand,  it  does  not  seem 
possible  to  establish  an  equation  for  F(s,t)  analogous  to  Eq.  (IV. 39), 


i.e.,  one  involving  5 only  via  > ( C ) • valid  for  all  £.  Even  if  such  a 


relation  were  to  exist,  it  would  be  necessary  to  determine  the  analytic 
form  of  XQU)  for  the  specific  process  under  consideration  in  order  to 
carry  out  the  required  integration  over  £. 

It  is  nevertheless  possible  to  conclude  that  for  t * t with  t( 


*c*  c 

defined  below,  the  mean  maximum  and  higher  moments  are  approximately  given  by 


00 

Z (t)  = rn  + n 1 d £ s"'1 


l-e-S^1 


(IV. 42) 


i.e., that  for  t z tc,  the  distribution  F(c,t)  in  Eq.  (II. 10)  can  be  replaced 
by  the  asymptotic  form  exp[->0U)t]  averaged  over  the  initial  distribution 
p(xQ).  The  time  t is  given  by  the  relation 


1 


*' = T>  "o)  = W 


(IV. 43) 


where  C0  is  the  smallest  value  of  t for  which  the  condition  (IV. 39a)  holds. 


We  assume  here  that  \QiO  is  a monotonical ly  non-increasing  function  of  5, 


which  would  appear  to  be  true  in  physically  interesting  cases.  Equation 
(IV. 42)  can  be  justified  as  follows.  If  t > tc,  then  for  £ « tQ,  it 
fo1  lows  from  the  definition  of  F(j;,t)  in  Eq.  (II. 7)  that  F(s,t)  « 1. 


Therefore,  F(c,t)  maxes  a negligible  contribution  to  the  integrand  £n_1[l  - FU.t)] 


of  Eq.  (II. 10).  The  detailed  form'  of  F(c,t)  in  this  region  of  integration 
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's  thus  unimportant.  It  ' : only  when  4 * that  the  ^orni  °f  Hf*»t) 
is  important.  With  CQ  defined  in  £qs.  (IV. *13),  (1V.38)  and  (IV. 39) 
then  allow  us  to  replace  F ( f; , t ) with  its  asymptotic  form.  Hence  if 
t i t it  is  ayain  sufficient  to  know  only  Q(0  = 1/T^(4)  to  obtain  the 
mean  maximum  and  higher  moments.  It  should  be  noted  that  for  fixed  t it£, 
the  first  moment  Z ^ ( t ) obtained  from  the  asymptotic  form  (IV. 4?)  is  expected 
to  be  more  accurate  than  the  higher  moments  Zn(t),  n >.  2.  This  is  because 
the  largest  contributions  to  Z ^ ( t ) come  from  small  values  of  4 where 
exp[ -X o (c ) t]  « 1 ,so  that  errors  introduced  through  use  of  the  asymptotic 
form  are  less  important  than  in  the  evaluation  of  the  higher  moments . 

It  should  also  be  observed  that  fur  large  t the  mean  maximum  Z^t)  is  a 
good  measure  of  tne  distribution  of  maxima.  It  can  readily  be  shown  via 
Eq.  (IV. 42)  that  the  dispersion  of  the  maximum  Z1 ( t ) tends  to  zero,  i.e. 

I 

This  is  'n  marked  contrast  to  the  mean  first  passage  time  ^(4)  which,  as 
shown  by  1 s dispersion  [Eq.  IV. 41a]  is  not  a good  measure  for  the  distri- 
bution of  first  passage  times  in  the  limit  as  £— ® • 

The  asymptotic  approximations  to  the  extremum  moments  Yn ( t ) for 
U tc  are  given  by  Eq.  (IV. 42)  with  r = 0. 

5.  THE  O.U.  EQUATION 

In  order  to  make  some  progress  beyond  the  above  formalism,  it 
is  necessary  to  consider  some  particular  Fokker-Pl anck  equation.  For  our 
first  example  we  turn  to  the  Ornstein-Uhlenbeck  (O.U.)  equation 


im 


z2(t)  - 


z,(t) 


V 


Zj(t)‘ 


1/2 


— 0 


(IV. 44) 
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r 


1 


at  f(x’^x0^  = ax 


I.x  f(x,t|xQ)  ♦ D f(x,t|x0)|  (V.1) 


where  the  general  coefficients  m^(x)  and  n^U)  of  the  F-P  equation 
(IV. 1)  now  take  the  form  rn^(x)  = - yx,  n^x)  * 2D  with  y,  D constant. 
After  defining  the  dimensionless  variables 


pt— ►! 


(p/2 D)1/2x-  y 


(V  .2) 
(V.3) 


we  can  rewrite  the  eigenfunction  equation  (IV. 7)  as 


d\(y) 

dy2 


dU.(y) 

- 2y  — * ♦ 2xkUk(y)  = 0 

dy 


(V.4) 


where  the  eigenvalues  xk  of  Eq.  (V.4)  are  dimensionless  ones  related  to 
these  of  Eq.  (IV. 7)  by  ,\k  - All  variables  and  functions  from  here 
on  are  understood  to  be  in  dimensionless  form, even  when  the  same  symbols 
are  used  for  them  as  in  earlier  sections.  The  weight  function  p(y) 


normalized  with  respect  to  y is 

.2 


L e-y 
\/7 


" /—  1 - erf  a 


(V.5) 


where  the  value  of  ^depends  on  the  boundary  conditions. 

We  now  consider  the  solution  for  each  type  of  boundary  condition 
separately. 

a)  Reflecting  Boundary  at  y = r,  Absorbing  Boundary  at  y = g. 

The  solutions  of  the  eigenfunction  equation  (V.4)  with  boundary 
conditions  (IV. 15)  and  (IV. 16)  are  combinations  of  confluent  hypergeometric 


functions : 


.(28) 


Uk(y)  = M 


/ Ak  I.  2 \ 2xkryH(1'7L,l,r  ) H(  7 ' 4»  h y2) 

T’"y  j i ^ y) 


(V.6) 
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The  eigenvalues  are  obtained  numerically  from  the  boundary  condition  (IV. 13) 
as  will  be  discussed  below.  Our  choice  of  normalization  makes  the  lowest 

eigenfunction  UQ (y ) — 1 for  any  finite  y as 

Equation  (V.6)  simplifies  in  a number  of  special  cases.  When  r—  - » 

( 28 1 

the  eigenfunctions  reduce  to' 


My)  = M 


,/-Xk  1 . 24 ' 4)  .Jl  ^ 3 2\ 

XT’’  2'*  )*  PkT^'2  2’2  ' 

rV  l) 


(V.7) 


where  r(y)  is  the  gamma  function.  For  the  special  case  r — - and  c 0 
the  eigenfunctions  are  the  odd  Hermite  polynomials,  Uk(y)  = H2k+1(y),  and 
the  eigenvalues  are  the  odd  positive  integers,  - 2k+l.  In  this  case 
the  full  probability  density  function  f (y ,t |yQ)  can  be  written  in  closed 


form  as 


f(y,i|y0>  ■ H'-e'"’)]  !«p  - 


- exp 


(y+y0e"T) 

(1-e  ) 


with  — < yQ  < 0.  The  cumulative  distribution  function  is 


-yoe 

F(  C*0,t  |yQ)  = erf  °_$x  "T/2 


(V  .8) 


(V.9) 


If  r _ and  t,  — °°  the  eigenfunctions  are  all  the  Hermite  polynomials, 

i 

U (y)  = H.(y),  and  the  eigenvalues  are  the  nonnegative  integers, 
k k 

. (29) 

The  probability  density  function  in  this  case  is 

..  i -2,.r'/2  f 


r -?T.r1/2  y'yoe 

f(y,x|y0)  = 1*0  - e )]  exp  ['^7^ 
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and  the  cumulative  distribution  function  F(c— ® | yQ)  is  of  course  unity. 

When  r s 0,  the  second  term  on  the  r.h.s.  of  Eq.  (V.6J  vanishes. 

If  in  this  case  °°,  the  eigenfunctions  are  the  even  Hermite  polynomials, 

U (y)  = H,.(y),  and  the  eigenvalues  are  the  even  nonnegative  integers, 
k 2 k 


c 2k.-  The  probability  density  function  is  as  in  (V.10)  with  -yQ  re- 
placing y , and  the  cumulative  distribution  function  is  again  unity, 
b)  Absorbing  boundaries  at  y = 

The  solutions  of  Eq.  (V.4)  with  boundary  conditions  (IV. 19)  are 
either  even  or  odd  functions  of  y: 


u2k(y)  ■ ^,y2)  , k - 0,1,2 

U2k+](y)  S yM(l  ‘ ?’y2)  ’ k = 0,1 ,2’"' 


(V.lla) 


(V.llb) 

Aaain  U (v)  - 1 as  c— • The  eigenvalues  are  obtained  from  the  conditions 


u2k<«>  * U2ktl(<>  = 0 


(V.12) 


The  odd  eigenfunctions  do  not  contribute  to  the  cumulative  distribution 
function  F(c,t)  with  0(y)  as  the  initial  distribution,  and  they  also  do  not 
contribute  to  FU,t|yo«0).  When  the  eigenfunctions  and  eigenvalues 

are  the  same  as  those  of  the  r — - case. 

Expressions  for  the  nvan  first  passage  time  T -j  (^ ) in  terms  of  the 
auxiliary  function  x(z)  of  Eq.  (IV. 26)  can  readily  be  evaluated  for  the 
equilibrium  weight  function  p(y)  of  Eq.  (V.5).  They  are. 


^(0  = rrri^TT  I dy  exp(y2)  (erf  y ‘ erf  r)C 


(V . 1 3) 


for  a single  absorbing  boundary,  and 

r 


T,U)  =y/r\  dy  exp  (y2)  erf2  y 
1 J0 


(V . 14) 
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for  two  absorbing  boundaries.  Note  that  Eq.  (V.14)  is 

identical  to  Eq.  ( V . 1 3 ) with  r = 0.  When  4 is  very  large,  asymptotic 
expansion  of  Eq.  (V.13)  y’elds^ 

2 

T 1 ( 4 ) - erf  r)  g-xP.iO..  . (V.15) 

c ) Resul ts  and  Discussion 

To  obtain  the  desired  results  for  more  general  boundary  conditions 
on  the  approach  of  the  Markov  properties  to  those  of  independent  random 
variables  it  is  necessary  to  calculate  numerically  the  eigenfunction  (y ) 
and  the  eigenvalues  x^.  These  calculations  were  performed  on  a Burroughs 
6700  computer  and  the  plotting  was  done  on  a CDC  3600  computer.  In  Tables 
I through  IV  and  Figures  3 through  7 we  present  a representative  sample 
of  the  extensive  calculations  that  were  carried  out.  The  complete  set  of 
calculations  is  available  upon  request. 


The  quantities  of  main  interest  are  tne  eigenvalues, since  all 
subsequent  results  depend  upon  them.  Since  p(y)  is  non-zero  for  all  finite 
y,  the  search  for  the  eigenvalues  reduces  to  finding  the  zeroes  of  the 
eigenfunctions  of  Sections  5a  and  5b.  A "bisecting  search"  method  was 
used  to  obtain  the  first  four  eigenvalues  as  a function  of  the  absorbing 
boundary  4,  for  various  values  of  the  reflecting  boundary  r and  for  symmetrical 
absorbing  boundaries.  The  results  for  some  representative  cares  are  given 
in  Tables  I to  III.  It  should  be  noted  that  the  ratio  A0(OAjU)  decreases 
rapidly  with  increasing  4.  For  4=1,  this  ratio  is  already  of  0(10'^) 
for  all  the  boundary  conditions  studied  by  us.  It  is  of  course  this  rapid 
separation  of  the  lowest  and  next  eigenvalue  with  increasing  4 which  causes 
the  asymptotic  results  to  hold  to  a very  good  approximation  for  finite 
values  of  the  boundary  value  4.  It  is  also  important  to  note  (see  Table  IV) 
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that  the  condition  (IV. 39a),  i t».  a = { p(y)dx  « 1,  holds  for  all  values 


of  C i 1 


Results  for  the  cumulative  probability  distribution  F(t,t) 


of  Eq.  (IV. 29)  are  presented  in  Table  IV  for  ref lect  ing  boundaries  at  r = 0, 
r = -1 , and  r * - for  various  times  x i . As  can  be  seen  from  the  equations 
and  discussions  inSection*  5a  and  5b,  the  cumulative  distribution  function 
F(£,t)  for  r = 0,  with  an  absorbing  boundary  at  y = is  identical  to 
F(£,t)  for  symmetrical  absorbing  boundariesat  y ■ * C.  The  entries  in  the 
first  row  of  Table  IV  thus  pertain  both  to  the  case  y=r=0,y=£  and 
to  the  case  y = - f, . The  results  presented  in  Table  IV  and  all  subsequent 
ones  discussed  here  were  obtained  by  using  the  first  four  eigenvalues 
Xq (O  through  X3(0- 

The  asymptotic  cumulative  distributions  F(C,t)-*-exp[-X0(f;)'l]> 

Eq.  HV.39),  are  also  presented  for  comparisons  in  Tahlp  IV.  Thp  "pxart." 
distributions,  obtained  with  the  use  of  our  eigenvalues  xQ  through  X^  and 
the  asymptotic  distributions  exp[->0(0T] .are  quite  close  in  value  (to  within 
15%  - 20%)  already  for  £ 3 1,  regardless  of  the  boundary  conditions.  For 
£ i 3 they  are  essentially  indistinguishable.3 

Results  for  the  mean  first  passage  time  and  mean  extrema  are  shown 
in  Figures  1-6.  Figures  1,  2,  and  3 display  the  natural  logarithm  of  the 
mean  first  passage  time  T-|(0  for  symmetrical  absorbing  boundaries  as 
discussed  inSection  *band  for  reflecting  boundaries  with  r = -1  and 


as  discussed  in^aV^he  solid  curves  are  the  exact  mean  first 


r = -® 


passage  times  obtained  from  numerical  integration  of  Eq.  (V.13).  It  ’s 


It  should  be  noted  that  the  values  of  F(1,t)  may  have  a non-negligible 
error  for  short  times  due  to  the  truncation  of  the  expansion  for  FU.t) 
after  four  terms. 


r — _z__ 

Iof  interest  to  note  that  the  mean  first  passage  times  obtained  from  the 
fifbt  four  terms  of  Eq.  (IV. 40)  agree  closely  with  these  exact  results. 


The  dashed  curves  are  the  independent  random  variable  approximation  1/^(0  e T,(r). 

For  £ £ 1 this  approximation  gives  results  which  are  extremely  close  to 
the  exact  T^C)  in  all  cases.  The  dotted  curves  represent  the  asymptotic 
expansion  ( V . 1 5 ) of  Eq.  (V.13).  It  is  rather  a poor  approximation  for  small 
values  of  £,  as  expected,  but  also  becomes  quite  accurate  for  £ > 1. 

We  have  also  evaluated  the  dispersion  of  the  mean  first  passage 
I 2W/2 

time,  (T^  - Tj  ) /T^.  It  is  approximately  unity  for  all  values  of  £ . 

This  extends  to  small  values  of  £ the  asymptotic  result  (IV. 41a)  on  the 
dispersion  of  the  mean  first  passage  time  for  the  general  Fokker-Planck 
equation.  This  is  an  interesting  and  disturbing  result.  From  it  one  learns 
that  the  mean  first  passage  time  T^(£)  is  not  a "sharp”  measure. and  one 
needs  really  to  calculate  the  full  first  passage  time  distribution,  Eq. 

(II. 9). 

The  mean  extremum  (T ) as  a function  of  time  are  shown  in  Fig.  4. 

By  extrenum.we  mean  here  the  greatest  excursion  of  the  random  variable  Y 
above  or  below  y = 0.  In  Figs.  5 and  6 we  display  the  mean  maxinum  Z.|  (x ) 
as  a function  of  time  for  processes  with  ref lecting  bounJaries  at  r = -1  and 
r = -«»,  where  by  maximum  we  refer  to  the  greatest  excursion  of  the  random 
variable  y above  y = 0.  The  dashed  curves  are  the  approximate  first  moments 
Z^(x)  computed  from  Eq.  (IV. 42).  The  solid  curves  are  the  exact  results, 
averaged  over  the  initial  conditions,  obtained  from  Eqs.  (II. 10)  and  (11.15) 
for  n = 1.  The  dotted  curves  are  the  differences  between  the  exact  and 
approximate  results.  It  will  be  noted  that  these  differences  are  very  small 
at  all  times  t and  essentially  zero  for  all  times  t > x in  agreement  with 
our  criterion  for  the  validity  of  (IV. 42).  From  the  definition  of  in 
terms  of  x0 ( £Q) , Eq.  (IV. 43),  one  can  see  from  Tables  I - 1 1 1 that  tc  is  of 
order  unity.  _n 


It  is  interesting  to  compare  this  value  of  tc  with  the  characteristic  time 
for  decay  of  correlations  in  the  unbounded  O.U.  process.  It  follows 
directly  from  Eq.  ( V . 1 ) in  terms  of  the  dimensionless  variables  of  Eqs. 

(V.2)  and  (V.3)  that  in  the  absence  of  boundaries 

on 

<y0  )y(0)>  E y„  j dy  y f(y,i  |y  ) = y 2 e‘T  . (V.16) 

characteristic  dimensionless  decay  time#i.e.  the  time  required  for 
the  correlations  to  decay  to  1/e  of  the  initial  value, is  thus  unity.  The 
values  of  xc  found  here  imply  the  physically  reasonable  result  that  mean 
extrema  and  maxima  of  the  O.U.  process  are  well  approximated  by  independent 
random  variable  results  at  times  when  the  correlations  of  the  unbounded 
O.U.  process  have  decayed  to  about  1/e  of  their  initial  value.  Analogous 
relations  for  more  general  Markov  processes  still  need  to  be  developed. 

We  have  also  calculated  the  dispersion  of  the  mean  maximum, 

2 1/2 

(^2~2i  ) /^i  • For  small  values  of  t it  is  of  0(1);  as  t + “ the  dispersion 

approaches  zero,  as  expected  from  Eq.  (IV. 44)..  In  contrast  to  the  relation 
between  the  mean  first  passage  time  and  the  first  passage  time  distribution, 
the  mean  maximum  is  thus  a useful  measure  of  the  distribution  of  maxima, 
particularly  at  large  times  t.  It  is  interesting  and  important  to  note 
that  one  can  therefore  make  a meaningful  calculation,  i.e.  one  with  a 
small  dispersion,  of  the  most  probable  maximum  excursion  of  the  space 
random  variable  for  any  time  interval  (0,i),but  that  the  calculation  of 
the  mean  time  required  for  the  space  variable  to  reach  an  extreme  value 
for  the  first  time  is  not  very  useful.owing  to  the  large  dispersion  of  the 
mean  first  passage  time. 
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Finally,  we  wish  to  point  out  that  the  units  in  this  discussion 
are  scaled  dimensionless  ones.  The  absolute  values  of  the  time  t and  the 
variable  x for  which  independent  random  variable  behavior  can  be  assumed 
depend  on  the  parameters  u and  D of  the  O.U.  equation,  according  to  Eqs. 

(V  .2)  and  (V.3). 

6.  HARMONIC  OSCILLATOR  DISSOCIATION 

A specific  example  to  which  the  developments  of  Section  4 can 
readily  be  applied  is  the  harmonic  oscillator  dissociation  model  of  Montroll 
and  Shuler  in  the  high  temperature  approximation. (®-10) 

Consider  an  ensemble  of  harmonic  oscillators  in  contact  with  a 
heat  bath  at  temperature  T;  the  fundamental  frequency  of  each  oscillator 
is  v.  Let 

0 h hv/kT  . (VI. 1) 

In  the  high  temperature  limit,  9 « 1,  one  can  approximate  the  discrete 

energy  levels  of  each  oscillator  by  a continuum  of  energies  denoted,  in 

units  of  hv,  by  the  dimensionless  variable  x.’  The  ground  state  energy  is 

(81 

x = 0;  it  is  assumed'  ' that  if  the  energy  of  an  oscillator  reaches  or 
exceeds  x = r,  the  oscillator  dissociates  irreversibly.  The  probability 
that  an  oscillator  has  not  dissociated  at  time  t then  obeys  a Fokker- 
Planck  equation  with  a reflecting  boundary  at  x = 0 and  an  absorbing 
boundary  at  x = 4 . The  coefficients  in  Eq.  (11.17)  are^’^ 

m1  (x)  = 1 - 0x  (VI. 2) 

n^x)  = 2x  . (VI. 3) 
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The  time  t is  expressed  in  units  of  1/k.jQ,  where  k^Q  is  the  collisional 
deactivation  rate  of  the  first  excited  state  to  the  ground  state. 

The  eigenfunctions  for  this  Fokker-Pl anck  equation  are  confluent 

, (9,28) 

hypergeometri<~  functions, 


Uk(x)  = M(-xk/0,l  ,0x),  k = 0,1,2,... 


(VI. 4) 


where  the  Ak  are  determined  by  the  absorbing  boundary  condition  Uk(0  = 0. 
When  £ - 00 , the  eigenfunctions  are  the  Laguerre  polynomials,  Uk(x)  = Lk(0x), 
with  eigenvalues  = k9.  The  weight  function  for  this  problem,  obtained 
from  Eq.  (IV. 5) , is 


p (x)  = 6 e" 


(VI. 5) 


The  mean  first  passage  time  to  x = t with  r = 0 (i.e.,  the  mean 


for  u ii>soc iaL iun  lu  unui  ) o given  by  Eq.  (IV. 25)  as 


VO 


(1  - e'*)2 


(VI. 6) 


In  the  limit  5 + °°,  it  then  follows  from  Eqs.  (IV. 39)  and  (IV. 41),  according 


to  which  T^(0  = [A0(0]  . that 

FU,t)  * exp[-02  £e'0r’t] 


(VI. 7) 


where  FU,t)  is  the  probability  that  T^O  < t.  The  results  in  Eqs.  (VI. 6) 

and  (VI. 7)  are  identical  to  those  obtained  from  the  rigorous  discrete 

(8) 

energy  level  treatment  of  the  Montroll -Shuler  model  in  the  high  temperature 
limit  as  N + 1 *£-*•»  [Eq.  VII. 6 of  reference  8],  It  should  be  noted  that 
these  results  were  obtained  here  with  much  less  effort  than  in  the  original 
papers, refs.  8 and  9. 
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Figure  7 shows  a comparison  of  the  exact  and  approximate  mean 
first  passage  times.  The  ordinate  is  in  units  of  OT^(C)  and  the  abscissa 
is  in  units  of  Of..  The  solid  curve  is  the  exact  mean  first  passage  time 
obtained  from  the  exact  expression  in  Eq.  (VI. 6),  the  dashed  curve  is 
0/A  (0  evaluated  numerically,  and  the  dotted  curve  is  the  asymptotic 
expansion  of  the  exact  result  of  Eq.  (VI. 6).  Beyond  Of,  = 3,  T^C)  = 1/A  (0 
to  within  11%.  For  9?,  = 3,  £ p(x)dx  = 0.05«1  ,so  that  the  criterion  o* 

Eq.  (IV. 39a)  for  the  validity  of  Eqs.  (IV. 39)  and  (IV. 41)  is  well  fulfilled. 
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APPENDIX  A 


We  illustrate  here  how  a particular  Fokker-Planck  equation,  namely  the 
Ornstein-Uhlenbeck  (O.U.)  process,  can  readily  be  obtained  by  starting 
with  a model  of  a Markov  process  defined  only  for  the  discrete  times  t i nAt. 

Our  starting  point  is  the  so-called  stationary  auto-regressive  process 
of  first  order 

X(tn+])  « a(At)X(tn)  ♦ A(tn+1)  (A.  1 ) 

where  X is  a random  variable  and  A is  a "noise"  term.  Equation  (A.l)  can 
be  rewritten  as 


X(tn  + At)  = a(At)X(tn)  + A( t n + At) 


(A.  2) 


We  now  go  over  to  continuous  time,  so  that  t becomes  the  continuous  variable  t. 

n 

This  permits  us  to  replace  the  noise  term  with 


t+At 

A(t  + At)  - f a(t)dt 
t 


(A.  3) 


Dividing  (A. 2)  [with  (A.3)]  by  At  and  taking  the  limit  At  -*•  0 yields 

(A. 4) 

Expansion  of  a(At)  in  a Taylor  series  about  At  = 0, 

(A.  5) 


dX(t)  _ lim 
dt  At+O 


-a-^  " ]-1  X(t)  + act) 

At 


i(  At)  = 1 + 

d ( At ) 


At  + 0[( At)  ] 

At:o 

and  substitution  of  (A. 5)  into  (A. 4)  then  gives,  to  order  At, 


= -/*X(  t ) + act) 

dt 


with 


AC  = . d sL&l 

d(  At) 


At=0 


(A. 6) 


(A. 7) 


Equation  (A. 6)  will  be  recognized  as  the  Langevin  equation.  If  we  now 
assume  that  &(t)  corresponds  to  delta-correlated  Gaussian  noise  with 


<a(t)>  * 0 

<dt)a{f)>  = 2os(t  - 1') 


(A. 8) 


where  D is  a constant,  then  it  is  easy  to  show  by  well  known  methods  that 
the  Langevin  equation  (A. 6)  integrates  to  the  Fokker-Planck  equation  (11.17) 
with  m^x)  =-/*x  and  m2(x)  = 2D.  The  stationary  auto-regressive  proces 
of  first  order,  Eg.  (A.l),  is  thus  equivalent  to  the  O.U.  process  [Eg.  (V.  )_1 
under  the  assumption  of  continuous  time  and  delta-correlated  Gaussian  noise. 


APPENDIX  B 

We  present  here  some  of  the  equations  used  in  the  numerical  computations 
for  the  O.U.  process. 

A useful  function  for  these  calculations  is 


vk(y)  -=  uk(y) 


(B.l) 


ax, 


where  the  right  hand  side  represents  partial  differentiation  of  the 
eigenfunctions  of  Section  V.a  and  V.b,  for  fixed  y.  Some  of  the  quantitir 
that  enter  the  numerical  computations  can  be  simply  expressed  with  th?  use 
of  this  auxiliary  function. 

Normalization  constants  N, 


Equation  (V  .4)  in  self-adjoint  form  is 


— o(y) 
dy  . 


dUk(y)  ] 
dy 


xkp(y)Uk(y)  = 0 


(B.2) 


Differentiation  of  this  equation  with  respect  to  Xk  yields  an  inhomogeneous 
differential  equation  for  Vk(y) : 


ay  [■<»> 


dVk(y) 

"dy 


+ xkp(y)Vk(y)  = -p(y)Uk(y) 


(B.3) 
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rl 


Multiplying  Eqs.  (B.2)  and  (B.3)  by  Vfc(y)  and  Uk(y)  respectively,  subtracting 
them,  integrating  ihe  diffe.w.ice  from  b to  C,  and  using  the  appropriate 
boundary  conditions,  yields 


\(0  ■ | 


[2o(OVk(OUk‘(0 


for  b = r 
for  b = -C 


(B.4) 


where  the  prime  indicates  differentiation  with  respect  to  y. 
Coefficients  At 

Integration  of  Eq.  (B.2)  between  the  boundaries  gives 


j dy  p(y)Uk(y) 


- °(b)u: 

V‘> 


(B.5) 


Substitution  of  Eqs.  (B.4)  and  (B.5)  into  Eq.  (V.3)  results  in 


Ak(ei  * -[»k(Ovk(E)T 


(B.6) 
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TABLE  I 


+ 

First  four  eigenvalues  for  the  O.U.  process  with  symmetrical  absorbing  boundaries  at  -f.. 
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1.9 
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2.9 
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3.1 

3.2 

3.3 

3.4 

3.5 
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•3.7 

3.8 

3.9 
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X (0 

Q 
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. 7984S0c*00 
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.25322Br4.o0 
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. 53 1 3 1 5c-0 1 
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.4954 1 Op-02 
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. 1 1 51 BOp-02 
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x1(0 
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. 132586F4-02 
. 96402 15  + 0 1 
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x2U) 
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. 78a 96 7 r * 0 r | 
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First  four  eigenvalues  for  the  O.U.  process  with  boundaries  at  r - -1 , y C 


f. 

x0(f.) 

x?U) 
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rirst  four  eigenvalues  for  the  O.U.  process  with  boundaries  at  r 


X0U) 

^(0 

c) 

L 

-1  .0 

.263720P  401 

.510">82c-40l 

. 75?666f4  0 

-0.) 

.23'*56'F40l 

.66557  IF 4OI 

.723666F40 

”0.8 

,219?38f401 

.66161  1 f'  4 0 1 

. 69r,6  6Br  4 0 

“0  j 7 

. 1 907005 < 01 

.9386905  4 u 1 

.6681  3 Jr  4 0 

-n.r, 

. 1 82 1 90r 4 0 l 

.6  1 6 2 2 0 r 4 0 1 

. 66 ] 635f 4O 

-0fS 

. 1666  36^.01 

. 396  7955  4O 1 

.615971 p*0 

-o.o 

. I5I5I6F4OI 

.3791965401 

.591 13sr+0 

-o;3 

. 1 376  26  F 4 0 1 

.3566195401 

.567 1 26r  4 0 

-0.? 

• 1 29  1 6 7 f 4 0 1 

,335666r401 

.56J935F40 

-o;i 

. 1 1 1677r40l 

.3173235401 

, 52 1 56 ?f  4 0 

oil 

.891066F400 

.2036755401 

.679266F40 

<v 

. 7897o7C+00 

.267/665401 

.659296540 

0.3 

. 69601  Of  4OO 

.?52008r40l 

.960161540 

0.6 

.609609v-*00 

.236653F  4O 1 

.92)7815+0 

0.5 

.530383r+00 

.2252765401 

.9062) O5+O 

o:& 

.658] 33f+00 

.21 2663546 1 

.307623F+0 

o;? 

. 3926621'-*  00 

.200813F+01 

. 37 ) 6 ) 5r  + 0 

0.8 

. 33367 1F400 

. 189716F401 

,3561 8?F+0] 

0:9 

. 260963F4  00 

.1793695401 

.3917185+01 

1 . o 

.236236F400 

. 1 69766F4  0 1 

,3280 1 9f+0 

1 . 1 

. 193179F400 

.1608555401 

,315081f+0] 

1.2 

. 15766/f,  JO 

.1526785401 

.30?398f40] 

1 .3 

• 1 ?6766f ♦ 00 

.1952055401 

.291667F40] 

* r H 

. 1 00636^  4 00 

. 1 3062'  401 

.nou  r835  + o 

1.5 

• 7873995-01 

• 1 3231 . 401 

.2708905+01 

li6 

.6066925-0 1 

.1268635401 

.261639F+0] 

1.7 

.659393F-G1 

. 1220685401 

• 25? 1 59r40 

U8 

, 36 1 35 1 p— C 1 

.1  17^955401 

.2959085+0] 

1.9 

•269672F-01 

. 1 1 6 ??6f  + 0 1 

.2383735+01 

2.0 

. 178817F-01 

.1111515+01 

.2320695+0] 

2.1 

. 126695F-01 

. 1085R65401 

. 226908c  + 0 ] 

2.? 

.362683F-02 

.1069895401 

,221963f+0] 

2.3 

. 580 ] 88T -02 

.1097975+01 

.217)605+0] 

?> 

. 381 882F-02 

. 1039725+01 

.213657F+01 

2^5 

.2658875-0? 

.1029565+01 

.2103635+0] 

2.6 

. 156860F-02 

.1016975401 

.207815540] 

2^7 

.9539665-03 

.1011665401 

.205762540] 

2.8 

,576027F-03 

. 1007535401 

.206168F+01 

2.9 

.3389305-03 

.1006835+01 

• 20?9 1 25+0 ] 

3.0 

. 19561?f-03 

.1003025+01 

.2019915+01 

3il 

. 1 10276F-03 

.1001865401 

.20)325540] 

3.2 

.6089965-06 

. 1001 1C5  + 01 

.2008585+01 

3.3 

• 329 1 8 1 F-06 

.1000665401 

.2005615+0] 

3> 

. 176  176F-06 

.1000365401 

•200331F+0] 

3.5 

.90221 1F-05 

.1000205+01 

. 200 1 98c  4 0 ] 

3.6 

.657567F-05 

. 10O01 1E  + 01 

.2001 155401 

3.7 

.227227^-05 

.1000065+01 

• 2000655  * 0 1 

3.3 

. 1 10500F-05 

. 1000035+0 1 

,20003654  0 

3:9 

.526216F-06 

.1000015+01 

.2000195+0] 

9lo 

.265693F-06 

.1000015  »oi 

.2000105+01 

-40- 

x3lt> 

.987000549) 
,9Sr>3t>0(  4.)  1 
.92  tf>0|i  -4  ')\ 
.89200  ‘>F*h  1 
.86279 U ♦ ')| 
.8336065  .M 
.805?^  + 9 1 
.777706V  + !)  1 
. 7509895  4 0 1 
, 725  0BHr  4>  n 1 
. 6 7 1 7 2 0 r 4f?  1 
. 652291^ 4 7 1 
. 62957 3r + M 
,60769  7r4  r;  l 

.5R6615r40  l 
.56632394^1 
,50681 7 F + 0 1 
. 528098 r 4 n 1 

.510150^4)1 
.69298 1 F + 0 1 
. 6 7658br4r,  l 
, O 6 0 9 5 n p 4 {;  i 
. 6 6 6 0 9 7 c 4 Q 1 
.631 999c  4 01 

• 6 1 8663f  4 o 1 
, 6 06085f  4 3 1 

.39A2O0F4G1 
.3831 99r4  0 1 

• 3723865  4 (w  1 
. 363320p4 ) 1 

• 35951 2F4  i)  1 
. 306666r4  0 1 
•3391225+01 
.332536r4il 

• 32667 1 F4  n 1 

• 321 520f  40 l 

.317060^401 

• 3 1 3266f  4 p 1 
,31 0095c  4 0 1 
, 307508F ♦ 0 1 
. 305996F  4 0 1 
» 3038975  + 0 1 
.302696F4  9 1 

,3017  br>F  ♦ 3 1 

• 30 1 1 6 5 f 4 0 1 
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TABLE  IV 

Exact  and  asymptotic  cumulative  distribution  functions  for  the  O.U.  process 


♦Note:  aU)  = r p(y)dy 


FIGURE  CAPTIONS 

FIGURE  1.  O.U.  process.  Mean  first  passage  time  to  symmetrical  absorbing 

boundaries  at  Exact  result,  . Independent  random  variable 

approximation, . Asymptotic  expansion  of  exact  result,  

FIGURE  2.  O.U.  process.  Mean  first  passage  time  tos  with  reflecting 

boundary  at  r * -1.  Exact  result,  . Independent  random  variable 

approximation, . Asymptotic  expansion  of  exact  result,  

FIGURE  3.  O.U.  process.  Mean  first  passage  time  to 5 with  reflecting 

boundary  at  r = -«°.  Exact  result,  . Independent  random  variable 

approximation, . Asymptotic  expansion  of  exact  result,  

FIGURE  4.  O.U.  process.  Mean  extremum.  Exact  result,  . 

Independent  random  variable  approximation, . Difference,  

FIGURE  5.  O.U.  process.  Mean  maximum  with  a reflecting  boundary  at  r = -1. 

Exact  result, -.  Independent  random  variable  approximation, . 

Difference,  

FIGURE  G.  O.U.  process.  Mean  maximum  with  a roflecting  boundary  at  r = -*>. 

Exact  result,  . Independent  random  variable  approximation, . 

Difference,  

FIGURE  7.  Harmonic  oscillator  dissociation.  Mean  first  passage  time  to  £. 

Exact  result,  . Independent  random  variable  approximation, . 

Asymptotic  expansion  of  exact  result,  
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f 


MEAN  EXTREMUM,  Y,  (t) 


MEAN  MAXIMUM,  Z,  (r) 


